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ABSTRACT 



New light curves of the gravitationally lensed double quasar Q0957+561 in the gr bands during 2008-2010 include densely sampled, 
sharp intrinsic fluctuations with unprecedentedly high signal-to-noise ratio. These relatively violent flux variations allow us to very 
accurately measure the g-band and r-band time delays between the two quasar images A and B. Using correlation functions, we obtain 
that the two time delays are inconsistent with each other at the 2<r level, with the r-band delay exceeding the 417-day delay in the g 
band by about 3 days. We also studied the long-term evolution of the delay-corrected flux ratio B/A from our homogeneous two-band 
monitoring with the Liverpool Robotic Telescope between 2005 and 2010*. This ratio B/A slightly increases in periods of violent 
activity, which seems to be correlated with the flux level in these periods. The presence of the previously reported dense cloud within 
the cD lensing galaxy, along the line of sight to the A image, could account for the observed time delay and flux ratio anomalies. 
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1. Introduction 



strong - black hole physics - galaxies: elliptical and lenticular, cD - quasars: individual 



The optical continuum variability of the gravitationally lensed 
double quasar Q0957+561 a t redshift z = 1.41 h as been widely 
studied since its discovery by Wal sh et al71(ll979l) . Several moni- 
toring campaigns focused on the determination of the time delay 
between the two quasar im a ges A and B (e.g ., Vanderriest et al. I 
Il989t IKundic et al.l Il997t [Serra-Ricart et al. 1999), where a 
major breakthrough occurred in Kundi cet al. I (1 1997b . who used 
Apache Point Observatory (APO) data. The 1.5-year monitor- 
ing programme with the APO 3.5 m telescope led to an accurate 
time delay AtgA = 417 ± 3 d (2<x confid ence interval; A lead- 
ing) in the g band. iKundic et al. I & 19971) also reported Af^A ~ 
420 d in the r band, which was consistent with the g-band delay 
measurement. A recent 2.5-year campaign with the Liverpool 2 
m robotic telescope (LRT) has confirmed the APO g-band de- 
lay, but it h as not allowed us to m easure a reliable time delay in 
the r band (Shalya pin et al. 1 12008. Paper I). Unfortunately, very 
accurate estimates of multiband delays between Q0957+561A 
and Q0957+561B remain elusive because of the absence of very 
prominent flux variations with signal-to-noise ratio S/N > 10, 
where S /N for a given fluctuation is defined as the ratio between 
its semiamplitude and mean flux error (see Paper I). The strong 
gravitational lensing scenario predicts the existence of an achro- 
matic delay (e.g., Schneider et al. 1992; Kochanek et al. 2004), 
while the possible detection of different delays in different opti- 
cal bands would provide extremely valuable information on the 
physical properties of the intervening medium. 



* Tables 1 and 2 corresponding to the Liverpool Robotic Telescope 
light curves are only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat 7J/A+ A/ 



Using the APO light curves for the two quasar images, 



Collier (2001) found that the r-band main fluctuations lag with 
respect to those in the g-band by 3.4 *j ^ d (l<x interval). This in- 
terband delay was interpreted as clear evidence for stratified re- 
processing within an accretion disc that is irradiated by a central 
high-energy source. The accretion disc would orbit the central 
supermassive black hole of the quasar. Interestingly, At rg ~ 4 d 
for the image B data alone, whereas At rg ~ 1 d for the image A 
data alone. These two estimates agreed within the lcr error bars, 
but the shortest delay from A data was thought to be underesti- 
mated as a result of the relatively poor sampling and variability 
behaviour dCollier 11200 IT) . The LRT follow-up of Q0957+561 A 
also led to interband delay estimates centred on 3^4.5 d (Paper 
I), seemingly supporting the four-day value for both images. 
We note that the presence of equal interband delays for the two 
quasar images is equivalent to the occurrence of equal delays 
between images in different bands. 

One can also obtain the delay-corrected flux ratio at time 
t: B/A = SB(t)/SA(t - AtBA), where Sa and Sb are fluxes 
of Q0957+561A and Q0957+561B, respectively. Although the 
strong gravitational lensing scenario produces achromatic and 
statio nary flux ratios of lens ed quasars (e.g., ISchneider et al. I 
1992; Ko chanek et al. I [20041) . actual scenarios are not so sim- 
ple. Chromatic flux ratios are usually related to differ ential 
extinction (e.g.. iFalco et al.l 1999; Eliasdottir et al. 20Qa) or 
differential microlensing (e.g., [Yoneharae t al. 2008). In addi- 
tion, time-variable flux ratios are likely due to differential mi- 
crolensing by stars in lensing galaxies (e.g.. Ilrwin et al. I IJ989: 
Para ficz et al. 1120061) . The light rays associated with the two im- 
ages of Q0957+561 pass through two separate regions within 
the central cluster cD galaxy at z — 0.36 acting as main gravi- 
tational lens (IStockt on 198ft lYounget al. I 1980; Ga rrett et all 
1992). Thus, while the optical continuum of the B image prob- 
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ably does not suffer significant dust extinction, the optical con- 
tinuum light of the A image is affected by a dense dusty cloud 
inside the cD galaxy (Goicoecheaet al. 2005a b). This differ- 
ential extinction produces a chromatic flux ratio, whose /?-band 
value w as basically constan t from 1987 through 2000 (e.g., see 
Fig. 3 of lOscoz et al. 112002b . The LRT observations also support 
the constancy of B/A over the 2000-2007 period in the g and r 
bands. Despite the fact that stars in the main lensing galaxy may 
induce variations in B/A (see above), the time-domain studies of 
Q0957+561 over two decades failed to detect these variations. 

We conducted a long-term photometric monitoring pro- 
gramme of Q0957+561 using the LRT at La Palma, Canary 
Islands. The observations are p art of the Liverpoo l Quas ar Lens 
Monitoring (LQLM) project dGoicoechea et al. I l2010h . Two- 
colour light curves of Q0957+561 during the first phase of this 
project (LQLM I; from January 2005 to July 2007) were pub- 
lished in Paper I. Here, in Sect. 2, we present new light curves in 
the g and r bands (LQLM II; from February 2008 to July 2010). 
These new light curves show densely sampled, sharp intrinsic 
fluctuations with S /N ~ 10, which are used in Sect. 3 to measure 
delays with unprecedented accuracy. In Sect. 4, we discuss the 
time-evolution of B/A over the 2000-2010 decade in the g and 
r bands. Our main conclusions are included in Sect. 5. In Sect. 
6, we briefly comment on some scenarios that could account for 
the observational results, as well as future prospects. 

2. Observations and data reduction 

All LQLM II optical frames of Q0957+561 were obtained with 
RATCam. This is a CCD camera with a 4'6 x 4!6 field of view, 
having a pixel scale of ~ 0727 (binning 2x2). To obtain a photo- 
metric signal-to-noise ratio of ~ 100 for the two quasar images 
for each observing night, we set the exposure times to 120 s per 
night in the g and r bands. Apart from the basic pre-processing 
tasks included in the LRT pipeline, we cleaned some cosmic rays 
and interpolated over bad pixels using the bad pixel mask. 

The pre-processed frames flow through our photometric 
pipelines to subsequent stages of processing (e.g., see the 
flowchart in Figure 1 of Goicoec hea et al. I 2010). At an initial 
stage, the crowded-field photometry pipeline produces the in- 
strumental fluxes of the quasar images. The frames that are of 
little or no interest were then removed from the initial data set. In 
Paper I we showed that quasar images with signal-to-noise ratio 
above 80 produce high-quality photometric results. Thus, only 
gr frames with signal-to-noise ratio > 80 over Q0957+561A 
are passed through the transformation pipeline. This pipeline 
transforms instrumental magnitudes into SDSS magnitudes, and 
the calibration-correction scheme is described in Appendix A of 
Paper I. 

For the long-term data in the g and r bands, we used a sophis- 
ticated transformation model incorporating zero-point, colour 
and inhomogeneity terms. Although this last term played an im- 
portant role when analysing earlier observations with the LRT 
(2005-2007; Paper I), the new data in the 2008-2010 period in- 
dicate that inhomogeneities have been weaker in the most recent 
years. Some improvements to the telescope in September 2007 
seem to have decreased inhomogeneities and typical seeing val- 
ues. To obtain g-SDSS magnitudes of a quasar image, we ini- 
tially considered an average colour {(g - r)sDss) m its colour 
correction (see Appendix A of Paper I). However, variations of 
(g - r )sDss ma Y introduce a non-negligible colour noise, which 
should be removed from the brightness record. The amplitude 
of the colour noises for the two images is ~ 10 mmag, and we 
eliminated these systematic noises in our 2008-2010 g-SDSS 
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Fig.l. Fluxes of Q0957+561A (top panel) and Q0957+561B 
(bottom panel) in the g and r bands of the SDSS photometric 
system. Circles denote the g-SDSS light curves and squares rep- 
resent the r-SDSS records. The whole optical data set is sepa- 
rated into two parts by a vertical dashed line: LQLM I (2005- 
2007) and LQLM II (2008-2010). The new effort in the 2008- 
2010 period doubles our contribution to the two-colour variabil- 
ity database of Q0957+561. We also highlight the LRT main 
fluctuations using grey rectangles (see main text). 



records. We also turned magnitudes into fluxes (in mJy) using 
SDSS conversion equations^. 

A large optical variability database, incorporating previ- 
ous LQLM I fluxes and the new LQLM II light curves of 
Q0957+561, is available in tabular format at the CDS0: Table 1- 
2 include 357 g-SDSS and 371 r-SDSS pairs of fluxes (S A ,S B ), 
respectively. Each of these tables contains the following in- 
formation. Column 1 lists the observing date (MJD-50000), 
Columns 2 and 3 indicate the flux and its error for the image 

A, and Columns 4 and 5 give the flux and its error for the image 

B. The LQLM optical light curves of Q0957+561 A (top panel) 
and Q0957+561B (bottom panel) are shown in Fig. Q] A verti- 
cal dashed line on 1 January 2008 separates the LQLM I and II 
periods. In the second monitoring period, there are 215 g-SDSS 
fluxes for each quasar image (circles), as well as 239 r-SDSS 
pairs of fluxes (squares). We achieve 1-1.3% photometric ac- 
curacy during the 2008-2010 period. Moreover, excluding the 
unavoidable seasonal gaps, the average separation between ad- 
jacent data is only three days. The new LRT light curves display 
four densely sampled, very prominent variations (see the two 
grey highlighted regions in Fig. [TJ. The two variations in S a are 
basically repeated in S b 14 months later, which means that these 



1 http://www.sdss.org/dr7/algorithms/fluxcal.html. 

2 Pipeline outputs, magnitudes and fluxes are also publicly available 
in the LQLM data-tools releases at http://grupos.unican.es/glendama/ 
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Fig. 2. LRT main fluctuations and some of their correlation functions. Top left panel: the fluxes of Q0957+561A (shifted forwards 
in time in 417 d; filled circles) and Q0957+561B (open circles) in the g band. Top right panel: the AB cross-correlation (solid line), 
and the average of the AA and BB autocorrelations (shifted by + 417 d; dashed line) in the g band for a — 6 d. Bottom left panel: the 
fluxes of Q0957+561A (shifted forwards in time in 417 d; filled squares) and Q0957+561B (open squares) in the r band. Bottom 
right panel: the AB cross-correlation (solid line), and the average of the AA and BB autocorrelations (shifted by + 417 d; dashed 
line) in the r band for a — 9 d. 



four fluctuations have an intrinsic origin. They also have S /N 
10, and can be considered as the LRT main fluctuations. 



3. Time delays from the LRT main fluctuations 

Although a g-band time delay between quasar images of about 
417 d is now firmly established, it is based on intrinsic flux vari- 
ations with 3 < S IN < 7 (see Paper I and references therein). 
Hence, the new g-band features within the grey rectangles in 
Fig.[T|(circles) represent a unique opportunity to measure the g- 
band delay with a very low uncertainty. These two well-sampled 
fluctuations with S IN - 13 are drawn together in the top left 
panel of Fig. [2] where S a (filled circles) is shifted 417 d for- 
wards in time. The time delay is determined by comparing the 
discrete cross-correlation function (DCF) to the discrete au- 
tocorrelation function (DAF). More properly, the delay corre- 
sponds to the minimum of the square difference between the 
DCF and the time -shifted DAF. This is the 6 2 technique (e.g., 
Serra-Ricart et al. 19990 relying on discrete correlation func- 
tions (Edelson & K rolik II 19881) . The 6 2 minimisation is a non- 
parametric method, which implies that one does not a priori as- 
sume a chosen model to relate the shapes of Sa and S b (see 
below). 

The differences Sa -(Sa) and S b~{S b) are the key pieces in 
the DCF, therefore we resampled both Sa and S b to obtain two 
curves with similar sampling (45 data points each), and to avoid 
biases between the averages (Sa) and (S b)- We evaluate the dis- 



crete correlation functions every day in two wide ranges of lags 
including correlation and anti-correlation peaks. The DCF and 
DAF were binned in 2a intervals centred at the lags, where a < 
10 d. For a — 3 d, both functions are very noisy, whereas for bin 
semisizes of 6 or 9 d, the discrete correlation functions in the g 
band have a smoother behaviour. In the top right panel of Fig. [2] 
using a - 6 d, we show the DAF shifted by + 417 d (dashed line) 
and the DCF (solid line). As expected, the two trends agree very 
well. We also followed a standard Monte Carlo approach to gen- 
erate 1000 synthetic data sets and determine time delay errors. 
In each synthetic light curve, the observed fluxes were modified 
by random Gaussian deviations that are consistent with the mea- 
sured uncertainties. We applied the 6 2 minimisation (see above) 
to each synthetic data set, and thus obtain 1000 delays for each 
value of a. Through the distributions of delays for bin semisizes 
of 6 and9 d, our final g-band measurement is AtsA = 416.5 + 1.0 
d (lcr interval). We also obtained the constraint: Af&4 < 418.5 d 
at the 99% confidence level. 

We repeated the procedure described in the two previous 
paragraphs, but using the r-band data in the bottom left panel 
of Fig. [2]instead of those in the g band. In this panel, the fluxes 
of the A image (filled squares) are shifted 417 d forwards in 
time. For a values of 6 or 9 d, the corresponding DCF and DAF 
are reasonably smooth. For example, the bottom right panel of 
Fig.|2]displays the DAF (shifted by + 417 d; dashed line) and the 
DCF (solid line) for a — 9 d. Surprisingly, the DAF should be 
shifted to the right by a few days to optimally match the DCF. 
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Fig. 3. Two-colour light curves of Q0957+561A. These curves 
correspond to slightly extended versions of the LRT main fluc- 
tuations of the A image: g band (circles) and r band (squares). 



To assess the significance of this extra delay (excess lag with 
respect to 417 d), we analysed in detail the delay distributions 
based on Monte Carlo simulations (see above). We find that the 
LRT r-band main fluctuations with S /N - 9.5 lead to AtBA = 
420.6 ± 1.9 d (lcr interval). Moreover, At BA is longer than 418.5 
d at 91-98% confidence levels, depending on the value of a. We 
can therefore state that chromaticity in AtBA is detected at about 
the 2cr level. 

This chromaticity in AtBA is supported by interband delays 
for the two quasar images. From the LRT main fluctuations in 
the g-band and r-band fluxes of Q0957+561B, we infer At rg > 
2.5 d at the 97% confidence level (a — 9 d). A very accurate, 
lcr delay At r „ =4+1 d was also obtained from these data 
dGil-Merino et al. 112012b . However, At rg < 2.5 d at the 98% con- 
fidence level (a - 9 d) for Q0957+561A data. This last con- 
straint is derived from densely sampled fluctuations in both op- 
tical bands (see Fig.O, which are slightly extended versions of 
the LRT main variations in the top panel of Fig. Q] It seems that 
previ ous claims of a four-day interband de lay for the two im- 
ages (Collier 2001; Shalyapin et al. 2008) were not accurate. 
Difficulties with relative ly low S/N val ues, monitoring gaps and 
other factors prevented ICollie r (2001) from separating a four- 
day delay for B from an one-day delay for A, and did not allow 
us to accurately determine At rg for A. Although we adopted an 
lcr interval of 4 + 2 d, some lcr lower limits in Table 2 of Paper 
I are equal or close to 1 d. 

The AB cross-correlation function is very sensitive to stan- 
dard microlensing variability, i.e., uncorrelated variations in the 
two quasar images. If the LRT main fluctuations (left panels of 
Fig.|2]i would be affected by microlensing, then their autocorrela- 
tion and cross-correlation functions would have different shapes, 
and the cr oss-correlation peaks w ould not reach a maximum 
value of 1 dGoicoechea et al.]|1998l) . However, these microlens- 
ing imprints are not seen in the right panels of Fig. [2] Although 
slow microle nsing was detected in light curves of several lensed 
quasars (e.g.. iGavnullina et al. 1 120051; iFohlmeister et al. I [2007; 



Shalvanin et al. 11200 9; Eulaers & Mag ain 11201 ll) . typical gradi- 
ents are too small to play a role in brightness records over rela- 
tively short time segments. For example, Hainl ine et al. I (1201 2|) 
used LQLM I data and more recent measurements from the 
United States Naval Observatory (USNO) to study the r-band 
flux ratio of Q0957+561. In this analysis conducted in paral- 
lel to ours, the authors report on the possible existence of a mi- 
crolensing gradient of 0.016 mag yr in the r band (see, how- 
ever, Section 4). Such a low gradient would produce an extrinsic 
variation of only 8 mmag over a six-month period, which is very 



much lower that the intrinsic signal that we find, and even lower 
than the noise level in our r-band data. Hence, the discussion 
throughout this paragraph indicates that standard microlensing is 
basically absent from the selected light curves and accordingly 
does not perturb the time delay estimates. 

Despite the robustness of non-parametric methods based on 
correlation functions, we also considered a x 1 (parametric) tech- 
nique to determine the time delay s between the two images (e.g., 
iKundic et ail 1 19971: lulfan et al. I f2006). This x 1 minimisatiorfl 
allows us to check the quality of the parametric model for relat- 
ing the shapes of S a and S b- The simplest model consists of a 
constant flux ratio, i.e., S s(t) = {B/A)S ^(f - AIba), where B/A is 
a constant. However, there is evidence for a B/A-S b correlation 
(see details and a discussion of this flux ratio behaviour in Sect. 
4-6), therefore we assumed an observationally motivated model: 
S B (t) = (B/A)S A (t - At^, B/A = (BIAUS B (t)IS BO f, where 
(B/A)o is the flux ratio at the reference flux S bo and B is the 
power-law index. Apart from the time delay AtBA, this scheme 
involves two additional free parameters (B/A)o and/?, which are 
used to link shapes. We do not know what the true way is to link 
S b to the correlated time evolution of S a- A power-law flux ratio 
is only one option among a variety of possible models, and there- 
fore, our^ 2 results should be taken with caution. To compare Sa 
and S b, we also used bins in A with semisize a. 

In the top panels of Fig. [4] we display AtBA-(B/A)o maps in- 
cluding our best solutions for a - 3 d (crosses) and their associ- 
ated 2cr contour lines. The g-band and r-band results are shown 
in the left and right panels, respectively. We obtained reduced 
chi-square values close to 1 (^- 2 /dof ~ 0.9), and two disjoint de- 
lay intervals around 417 d (g band) and 424 d (r band). These 
g-band and r-band delay intervals are separated by 4 d, whereas 
the difference between the best solutions is 7 d. Thus we find that 
the chromaticity in AtBA is more pronounced than that from the 
6 2 method (see above). Although the results for a - 3 d support 
a significant chromaticity of the delay, other values of a produce 
overlapping delay intervals. For example, if we take a — 9 d, the 
best solutions are characterised by ;f 2 /dof ~ 1.2-1.5. These ap- 
pear in the bottom panels of Fig. [4] (crosses). We also show the 
2<x contour lines around 418 d (g band; bottom left panel) and 
420 d (r band; bottom right panel). For a — 9 d, both delay inter- 
vals overlap with each other, and the parametric technique does 
not separate the two delays in the two optical bands. However, 
the production of "excessive chromaticity" or achromaticity is 
not surprising, since the degeneracy between the two shape pa- 
rameters and the delay likely prevents accurate/reliable^ 2 -based 
delay measurements. 

Hereafter, we consider the self-consistent solution for the de- 
lays: (a) At B A - 417 d in the g band, (b) At B A - 420 d in the r 
band, (c) At rg = 4 d for the B image, and (d) At rg — 1 d for the 
A image. Our solution agrees with the discussion in the previ- 
ous paragraphs of this section. We also note that several stud- 
ies of AtsA in the red part of the optical spectrum favoured de- 
lays above 417 d that are only margi nally consistent with the 
APO 2cr interval in th e g band (e.g.. ISerra-Ricart et al.ll!999l 
lOvaldsen et al. Il2003al) . This discrepancy was not originally as- 
sociated with a chromatic delay between images, but with less 
quasar variability and greater contamination from the lensing 
galaxy in red filters, the existence of multiple achromatic delays 
in long-term light curves, etc. 



3 We indeed minimised ,^ 2 /dof, with 'dof being the degrees of free- 
dom 
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Fig. 4. Best solutions and 2cr contour lines from the x 2 technique. Top left panel: a — 3 d (g band). Top right panel: a — 3 d (r band). 
Bottom left panel: a — 9 d (g band). Bottom right panel: a — 9 d (r band). 
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4. Two-COlour flux ratio over the 2000-201 decade Table 3. Flux ratio in the g band. 

The spectral behaviour and the long-term evolution of the 
delay-corrected flux ratio B/A has attracted incre asing atten- 
tion in the first decade of this century (e.g.,[ Refsdal et al. 2000; 
Oscozetal.1l2002t lOvaldsen et"al~T l2003bt iGoicoechea et al. I 
2005a b). The 1999-2000 Hubble Space Telescope (HST) spec- 
tra of Q0957+561 indic ated the chromaticity of B /A at optical 
continuum wavelengths (Goicoec hea et al.1l2005al) . At the aver- 
ag e wavelengths of the g an d r bands, the HST data in Fig. 1 
of IGoicoechea et al. I (l2005al) lead to lcr intervals B/A - 1.10 ± 
0.01 (g-band) and B/A = 1.04 + 0.02 (r band). It can be also 
demonstrated that the C in] (A 1909) and Mgn (A219S) emission 
lines only slightly influence the estimation of the optical contin- 
uum flux ratio from the g and r broad filters, introducing a small 
bias of (-B/A) cont /(B/A) cont+ ii ne ~ 1.01. Apart from spectral anal- 
yses, time domain studies suggested the c onstancy of the flu x 
ratio between 1987 and 2000 (R band; e.g. JOscoz et al.l l2002). 
and then from 2000 to 2007 (g and r bands; Paper I). Here, the 
new LRT data allow us to discuss the gr flux ratio in the 2007- 
2010 period, as well as to compare it with the HST two-colour 
ratio in 2000. 



" See main text. 

b Bin semisize. 

c dof = degrees of freedom. 

d 2a confidence interval when x 1 /dof < 1.2. 



To evaluate the flux ratio in the g band, one should compare 
the light curve of the B image and the fluxes of the A image 
shifted by + 417 d. For example, the Q0957+561B fluxes be- 
tween day 3649 and day 3894 have a very short counterpart in 
the original light curve of Q0957+561A before day 3477 (see 
the circles in both panels of Fig. [TJ. The counterpart only con- 
sists of three data points, and we did not calculate the g-band 
flux ratio over days 3649-3894. We also emphasize the absence 



of a counterpart in the r-band fluxes of A for this time segment 
of B. The first useful time segment of B covers days 4010^4249, 
and it is labelled TS 1 . The fluxes of B in TS 1 have a relatively 
long overlap with fluxes of A shifted by + 417 d (see the mid- 
dle and bottom left panels of Fig. 7 in Paper I). We removed 
two data points from the overlapping record of A because these 
fluxes are affected by atmospheric and/or instrumental problems 
(see Paper I for details). The other useful time segments of B are: 
TS2 (from day 4498 to day 4627), TS3 (from day 4777 to day 
4992) and TS4 (from day 5107 to day 5373). 

We used a^ 2 minimization to find the flux ratio for each seg- 
ment. In general, the shifted epochs of A do not coincide with 
the epochs of B, therefore we introduce bins in A around the 
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Fig. 5. AB comparisons in the g band. We show the overlapping periods between the light curves of A (filled circles) and B (open 
circles), where the fluxes of A are shifted by + 417 d and are properly amplified (see main text). Top left panel: the overlap period 
for the first time segment of B, i.e., TS1 (days 401CM-249). The A signal is amplified by 1.086. Top right panel: the overlap period 
for TS2 (days 4498-4627). The A signal is amplified by 1.14. Bottom left panel: the overlap period for TS3 (days 4777-4992). The 
A signal is amplified by 1 .086. Bottom right panel: the overlap period for TS4 (days 5 107-5373). The A signal is amplified by 1 . 14. 
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Fig. 6. AB comparisons in the r band. We used the original epochs and fluxes of B (open squares), while the fluxes of A are shifted 
by + 420 d and properly amplified (filled squares). Together with these LRT data, we also show the more poorly sampled USNO 
fluxes of A (filled triangles) and B (open triangles). Top and bottom left panels: comparisons in two periods of normal activity. The 
A signal is amplified by 1.019 (see main text). Top andbottom right panels: comparisons in two episodes of violent activity. The A 
signal is amplified by 1.057 (see main text). 
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Table 4. Flux ratio in the r band. 



Time segment^ 


c£(d) 


Best fit 


;r7dof 


B/A± 


TS1 


6 


1.023 


50.5/46 


1.019-1.027 




9 


1.023 


59.7/50 


1.019-1.026 


TS2 


6 


1.054 


52.0/46 


1.050-1.058 




9 


1.056 


80.5/51 


- 


TS3 


6 


1.016 


39.3/48 


1.012-1.019 




9 


1.016 


34.1/50 


1.012-1.019 


TS4 


6 


1.060 


160.5/54 


- 




9 


1.060 


169.0/57 


- 



" See main text. 

b Bin semisize. 

c dof = degrees of freedom. 

d 2tr confidence interval when ^ 2 /dof < 1.2. 



epochs of B (e.g.. lUllan et al. 1120061) . These bins have a semi- 
size a of 6 or 9 d. In Table [3] we give the best solutions and 
their reduced chi-square values, as well as some 2<x intervals 
for B/A. A constant flux ratio B/A - 1.086 works on both TS1 
and TS3, and it agrees well with the corresponding HST ratio 
in 2000 (1.086 x 1.01 ~ 1.10; see above). The AB comparisons 
for the time segments TS1 and TS3 are shown in the top and 
bottom left panels of Fig. |5J where we amplified the A signal by 
a factor of 1.086 (filled circles). Curiously enough, the simplest 
scenario (constant B/A) does work on TS2 and TS4, since most 
of the best solutions are associated with reduced chi-square val- 
ues ranging from 1 .5 to 2.8. These best solutions B/A ~ 1. 14 also 
differ from the ratio for TS1, TS3 and the HST observing dates. 
Taking an amplification of 1 . 14 for the signal A in TS2 and TS4, 
both A and B signals are compared to each other in the top and 
bottom right panels of Fig. H](A = filled circles and B = open 
circles). The flux ratio seems to reach "anomalous" values dur- 
ing episodes of violent activity (involving flux gradients > 0.1 
mJy/100 d; TS2 and TS4), while it remains lower and basically 
constant in other periods with flux gradients < 0.05 mJy/100 d. 
In the bottom right panel of Fig. [5] we also find evidence of a cor- 
relation between flux ratio and level of flux in TS4, i.e., B/A < 
1.14 at S B = 0.42 mJy, B/A ~ 1.14 at S B ~ 0.5 mJy and B/A > 
1.14 at S b = 0.59 mJy. This kind of correlation is not so evident 
in TS2. 

In Table |4] we present our results in the r band. For the two 
periods of normal activity in TS1 and TS3, a constant flux ratio 
B/A - 1.019 can account for the flux gap between the light curve 
of B and the record of A shifted by + 420 d (see Table [4] and the 
left panels of Fig.©. Although the adopted solution B/A - 1.019 
is only marginally consistent with the analyses in both periods, 
this should not cause suspicion of a possible decrease of B/A. 
The best-fit ^ 2 /dof values for TS3 are 0.7-0.8, so the formal 2cr 
intervals for this time segment only contain reduced chi-square 
values equal to or less than 0.9. Consequently, there are solutions 
B/A ~ 1.021-1.022 with^ 2 /dof ~ 1, which suggest that the flux 
ratio uncertainty for TS3 in Table His underestimated. Taking 
the small perturbation by the Mg n (/12798) line into account, we 
obtain a continuum flux ratio of B/A ~ 1.019 x 1.01 ~ 1.03 
at the average wavelength of the r band (see above). This ratio 
agrees with the HST determination of B/A at the same wave- 
length. For the two periods of violent activity (TS2 and TS4), a 
constant flux ratio does not convincingly explain the LRT obser- 
vations. Moreover, the best solutions B/A ~ 1.057 do not agree 
with those for TS 1 and TS3. We display the AB comparisons for 
TS2 and TS4 in the right panels of Fig. [6] where the open squares 
trace the light curve of B, and the delay-corrected and amplified 



fluxes of A are represented by filled squares. Once again, we find 
evidence of a B/A - S b correlation in the bottom right panel of 

Fig. |6l but this time in the r band. 

As we commented in Section 3, Hainli ne et al. I (120 12) found 
a slow gradient in the r-band flux ratio (in magnitudes) of 
Q0957+561. They used published LRT magnitudes together 
with new USNO data covering a general period similar to ours. 
However, we do not detect any long timescale drift in our anal- 
ysis with only LRT data, and this discrepancy needs more at- 
tention. First. lHainiineetal.l d2012l) adopted a LRT-USNO pho- 
tometric offset of 14.455 mag, which seems to be biased in + 
0.025 mag when comparing LRT and USNO magnitudes at sim- 
ilar epochs. In their Fig. 2, the authors derive the flux ratio in 
the second time-segment (days 4100^-200; it corresponds to our 
TS2) from a few differences Am,t(LRT) - Awb(USNO). Because 
the USNO fluxes are likely underestimated in 0.025 mag, these 
differences should be enlarged until reaching the values for the 
fourth time-segment (days 4750-5000; TS4 in our framework). 
Second, the flux ratio in the third time-segment (days 4500- 
4600; TS3 in our framework) is inferred from only a few USNO 
magnitude differences that include some outlier. In Fig. [6] we 
compare the LRT and USNO fluxes using an unbiased LRT- 
USNO photometric offset of 14.43 mag, as well as turning HJD 
into MJD and magnitudes into mJy. The open triangles describe 
the USNO light curve of B, and the delay-corrected and ampli- 
fied USNO fluxes of A are displayed as filled triangles. The time 
delay and the amplifications are those obtained from the LRT 
data (see above). In general, the LRT and USNO data agree very 
well. However, there is a clear outlier in the USNO record of B 
for TS3. This noticeable deviation from the general trend means 
that Sb is overestimated, and therefore, its associated magni- 
tude should be increased by a certain amount. In Hainline et al.'s 
scheme, this would lead to a lower Am^ - Am^ value in the third 
time-segment, so the new cloud of magnitude differences would 
more closely resemble the cloud in the first time-segment (TS 1 
in our framework). Therefore, both the LRT and USNO data sets 
are consistent with an oscillating behaviour of B/A. 



5. Conclusions 

Our main conclusions are: 

1. New LRT light curves of Q0957+561A and Q0957+561B in 
the gr bands during 2008-2010 show well-sampled, sharp 
intrinsic fluctuations with S/N ~ 10. These extraordinary 
features allowed us to very accurately determine the g-band 
and r-band time delays between both quasar images. The two 
time delays are inconsistent with each other at the 2<x level. 
More specifically, while we obtain a delay near to 417 d in 
the g band, there is an extra delay of about three days in the 
r band. This extra delay cannot be attributed to low S/N val- 
ues, contamination from the main lensing galaxy, or similar 
artifacts. 

2. From the LRT two-colour photometry of Q0957+561 dur- 
ing 2005-2010, we inferred the g-band and r-band delay- 
corrected flux ratio in four different time segments. The flux 
ratio B/A has an oscillating behaviour in both optical bands, 
reaching higher values in the two segments of violent activ- 
ity and remaining lower in the other two periods of normal 
activity. These normal activity periods are characterised by 
g-band flux gradients < 0.05 mJy/100 d, and B/A in each 
band does not vary from period to period or within a given 
period. The normal g-band and r-band ratios are also consis- 
tent with the HST ratios in 2000 at the average wavelengths 
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of the g and r bands. For the two episodes of violent activity 
(showing g-band flux gradients > 0.1 mJy/100 d), the flux 
ratio in each band is similar in both segments, but it seems to 
be correlated with the intra-segment flux level. 



6. Discussion and future work 

The optical continuum of Q0957+561A is plausibly affected by 
a dense dusty clou d inside the cD lensing galaxy at z — 0.36 
dGoicoechea et al. 1l2005a,b). Because the light propagation time 
in the intervening medium is expected to increase with decreas- 
ing wavelength (chromatic dispersion; e.g jBorn & Wo lf 1999), 
the presence of this substructure could be responsible for a three- 
day lag between g-band and r-band signals, and thus explain 
the observed chromaticity of the delay between images. A de- 
tailed discussion on the composition and size of the cloud along 
the line of sight to the A image is beyond the scope of this 
paper. If this scenario turns out to be true, it would be nec- 
essary to estimate the proper delay to obtain a refined delay- 
based determination of th e Hubble constant Hq (e.g.. iJacksonl 
120071; iFadelv et al. ll20Toh . In addition, future multi-wavelength 
(optical) monitoring campaigns of other gravitationally lensed 
quasars may also lead to unexpected delays caused by substruc- 
tures in non-local lensing galaxies, and thus, to improved esti- 
mates of lensing mass distributions and Hq, 

There is at least one crude interpretation for the flux ratio 
anomaly during violent episodes in Q0957+561. The violent ac- 
tivity may be related to a strong outflow, inducing a significant 
polarizati on degree in the otherwise weakly polarised UV emis- 
sion (e.g.. iBegelman & Sikoralll987b lBeloborodovlll998l) . For 
the A image, this polarised radiation would pass through a dust- 
rich region with alligned elongated dust grains, suffering from a 
higher extinction than that observed in periods of normal activity 
(dichroism; e.g.. lBorn & Wolflll999l) . The induced polarisation 
degree could increase with increasing activity of the central en- 
gine (flux level), so that more extinction would be observed for 
higher fluxes. Future polarimetric data in both normal and vio- 
lent periods will be used to check this interpretation. 

A microlensing scenario is difficult to reconcile with ob- 
servations of Q0957+561 for the last 25 years. The analysis 
of the LRT light curves indicates the absence of uncorrelated 
variations in the two quasar images (standard microlensing), 
and only a slight B/A increase occurs for the sharpest intrin- 
sic events. To account for this flux ratio anomaly, one might in- 
voke the possible existence of radial expansions of the accre- 
tion disc during violent episodes (a model of an accretion disc 
with a time-varying size has a lso recently been proposed by 
Blackbur ne & Kochanek 11201 Oh . The expanded sources would 
cover larger regions of the microlensing magnification pattern 
for the B image, and produce slight extra magnifications of that 
image. Although such an exotic microlensing seems t o work, the 
"exce ssive constancy" of B/A over ~ 25 years ("e.g. JOscoz et al. I 
2002, and this paper) calls this scenario into question. We think 
the next logical step should be to accurately study light curves 
covering more than 5-6 years. New 1999-2005 IAC-80 data in 
the R band3 together with the 2005-2010 LRT and 2008-2010 
USNO data in the r band will make up a 10 -year variability 
database, whereas additional old /?-band fluxes dOvaldsen et al. I 



l2003atlSerra-Ricart et al. Ill999l) and 201 1-2012 frames in the r 
band may contribute to a 20-year baseline. 
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(UC) collaboration. These frames will be fully reduced in a near future. 



